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ABSTRA.GT 


Demand for light weight aircraft structures results in the use of 
as small a safety margin as is practical. As a consequence of the 
small safety margin and other uncertainties, cracks or partial damages 
are likely to occur before the economical life of the aircraft is 
expended. Fatigue is one of the principal causes for the cracks. 

Fatigue loading and fatigue crack growth also contain uncertainties. 

The susceptibility of the aircraft structure to crack or partial 
damage during the useful life of the structure imposes the requirement 
that the structure should he capable of supporting the service loads 
with these cracks. Furthermore, it must be possible to detect these 
cracks before they extend to critical si^es and cause .catastrophic 
failure of the structure. Therefore, any fail safe design that can 
achieve this objective needs a knowledge of the probability of the 
presence of a crack of a certain length at a given location after 
certain number of flight hours. A stochastic model has been developed 
to describe such a probability for fatigue process by assuming a varying 
hazard rate. This stochastic model can be used to obtain the desired 
probability of a crack of certain length at a given location after certain 
number of cycles or time. 

Quantitative estimation of the developed model has also been 
dLscussed. Application of the model to develop a procedure for reliability 
based cost-effective fail-safe structural design has been discussed. 

This design procedure includes the reliability improvement due to 
inspection and repair. Methods of obtaining optimum inspection and 
maintenance schemes have also been discussed. 



Alt&rnate methods of fatigue reliability improvement by cold 
working processes have been discussed. The associated stress corrosion 
problem has been studied. Application of statistical decision theory 
to select suitable test options and safety factors subject to a reli- 
ability constraint have also been investigated. 

Most of the investigations under this project have either been 
published in journals and conference proceedings or pending publi- 


cation. 



INTRODUCTION 


It is now generally accepted that all structural materials are not 

"flaw- free*’. Sometiraea, a maximmn acceptable flaw-siae can be specified 

as a part of the structural specifications. Thus, an Initial flaw size 

(a ) and the associated probability distribution characterize the 
o 

strTicture. Due to fatigue loading, these initial "micro- sized” flaws 
groxf to detectable sizes. Time or number of cycles required for this 
groxirth to detectable size of crack length is often called the "crack 
initiation time.” In many cases, this grox?th time amounts to a 
significant function of the total fatigue life of the structure. Due 
to further fatigue loading, crack sizes increase until they attain critical 
sizes. These critical sizes depend upon the critical stress intensity 
factors and the external loads. Thus, the probability distribution for 
crack sizes is changing continuously x-Tith time or number of cycles at 
all locations of the structxire. Therefore, the probabilistic description 
of fatigue process can be e^pJ^essed as the probability that for time t ^ T, 

the crack size a ^ A. This is a shochastic process. 

X "8 

In most of the reported x^orks , the varying crack lengths "a” 
associated xTith the fatigue process are ignored. In these works, 
stochastic process is not considered. The entire fatigue process is 
described by a single random variable ”t" xdiich is the time for fatigue 
failure. The quantitative description consists of the probability that for 
time t < T fatigue failure xijill take place. Because of the simplicity of the 


model, probability distributions sucb as the Weibull Distribution have 
been used to describe the time for fatigue failure. 

The use of such a description that relies on a single random variable 
is very limited because such a model neglects many important aspects of 
fatigue process. For example, one question that needs an answer is as 
follows. What is the length of the crack that corresponds to the defined 
failure time? Is this length the initiation length or critical crack 
length or some arbitrarily chosen length? Initiation length can vary 
depending on the available non-destructive inspection capability. 
Furthermore, such a model does not provide any information for optimising 
repair threshold crack length,- crack arresting devices, W.D.I. capabilities 
and different loading process. Another arguement used by the users of a 
single random variable is to assume that the effect of varying crack 
length is negligable and a stochastic process is not needed. In order 
to verify if such a statement could be true fatigue data from specific 
fleet of aircraft are analyzed as a first step of the investigation. As 
explained in ^ater sections, these investigations demonstrated that a 
stochastic model is necessary to describe fatigue of structure . 

Further investigations during the project period are described as follov7s: 

a) Development of a simple stochastic model for fatigue by using 

9-10 

the concept of a varying hazard rate and a birth process , 

b) Quantitative estimation of the parameters of the stochastic 

9 

model by using fracture mechanics considerations^, 

c) Quantitative estimation of the parameters of the stochastic 

' , , ^ 11,12 

model from available data . 
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As a by-product o£ the above investigations^ an improved mathematical 
technique has been developed. This technique and its application can be 
described as follows; 


i) An improved numerical technique of multiple integration with respect 

17 

to one independent variable . 

j) Application of new technique of integration to develop a procedure 

18 

for the study of some random vibration problems . 
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Analysis of Fatigue Failure Data 


In order to investigate 1£ fatigue provess can be described by a single 
random variable "t" that denotes time for fatigue, fatigue data from two 
specific fleets were analyzed. A typical inspection record contained the 
following information. 

1. Identification number of the airplane 

2. Mumber of flight hours completed before the inspection 

3. Inspection date 

4. Number of reinspection (s) 

5. The command 

6. The base 

7. Facility of inspection 

8. Crack location by numbers of the critical regions as has been 
previously identified 

9. Number of such cracks in a given region 

10. Direction of crack growth 

11. Crack length 

12. Information as to ^diether the crack has been repaired 

A two parameter weibull model was hypothesized for the fatigue failure time ”t", 

F^ir) = t- [-(*/a)V ('} 

In the equation, a and 0 are shape and scale parameters respectively. These 
parameters were estimated from the data by using the method of maximum likelihood 


The chi-square and Kolmogrov tests were used to verify the goodness of 
fit of the estimated parameters. The following conclusions were reached, 

(1) For a given critical location or a selected group of critical 
locations, no acceptable Weibull distribution was obtained 
unless the data were censored in some way. In general, 
censoring of both high-level outliers and low level outliers 
were needed. low-level outliers refer to those fatigue failure 
times that lie below a selected failure time for purposes of 
censoring. Similarly, high-level outlier refers to those 


failure times that lie above a time corresponding to high« 
censoring level. Use of a low level outlier censor was not con'- 


servative.' Any model derived by the use of low level censoring 
can result in serious errors in decisions concerning design and 
maintenance. Similarly, models derived by the use of high level 
censoring caii result in increased weight and cost. 

(2) A three parameter weibull distribution or a log-normal distribution 
did not improve the results. 

(3) However, when the observed failure times at a given location 
were reduced by regression techniques to correspond to the time 
for initiation of crach of a given length acceptability of the 
two- parameter weibull model improved in many cases. Probability 
distribution 'Was different for different;-. crack lengths. 

Necessity for a stochastic model was evident. 




DavalopmEnt of A Stochastic Modal for Fatigue 


It is assumed that a single crack is present in a fatigue critical 
region. Multiple cracks can be treated by order statistics or other 
procedures. Thenj the variation of crack length with time is quali- 
tatively of the type shoxm in Figure 1, This consists of a continuous 
variation of crack length with continuous variation of time or number 
of cycles. The corresponding model for the stochastic process for 
fatigue crack siaes involves the consideration of continuous state 
space of crack lengths and continuous time. It is difficult to develop 
such a model. The development of the model is simplified by considering 
the state space of crack length to be discrete as shown in Figure 1. 
Accuracy can be increased by decreasing the magnitude of of 

discrete crack length increments. This process of considering the 
state space of crack length can also accomodate consideration of crack 
initiation i.e, probability of a crack of length a^ initiating at time t 
less than or equal to as shorn in Figure 2. Even though the crack 
lengths are assumed to increase in discrete steps the mean crack growth 
rate can vary continuously as a function of time. Because the resulting 
process is nonstationary, the probability that a crack of length k(i^0.> 
i.-.e. , -fe. times Mr" , is present at a time t <t^ depends on the initial 
value of time t . This is denoted by P(k, t , t, ) . 

By considering the different ways in which the event of the develop 
ment of a crack of length can occur in time interval t^ to t + 

the following equation can be written 
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^ p:,(k,*c,i)p , t+At) (zj 

i S.O 


by assuming ordetliness of crack growth i.e., 


and 


^ o 

C'3) 


LaV- S( 

= 0 
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where 

hi 

P(C,tc^ t) 

(s) 

and 

Sjt) =£ P(^,k‘>,i) -- 

1 — to, ) — P 

(&) 


lifl 

It: can be shown that the following differential equations are for k > 1 


In this equation j E [a(t)] is the mean crack growth rate at t. For 
k = 1, the equation (7) takes the following special form. 

^ P(htro,t) (g) 

where f (t) is the probability density for crack initiation. These 

G 

equations can be solved by methods similar to those discussed in Reference 
C^0« However, F(k, t) can be obtained only if E fa(t)] , f (t), 
and the probability distribution for initial flaw sizes are 
known. This procedure will be discussed in a separate note. The method 
of obtaining E[a(t)] is discussed in the next section. 

Mean Crack Growth Rate 

Knowledge of the mean crack growth rate is essential to estimate 

■ 21 

the crack length at a given time. According to Forman the rate of 
crack growth is given by 
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cJU^ 


■V) 


c 


, (aK) 






where G^ and n are material constants, illC is the range of stress intensity 
factors, is the critical-stress intensity factor, r is the ratio 
of minimtam stress intensity factor to the maximum stress intensity factor 
'a' is the half crack length and 'N' is the number of cycles. For a 
stiffened panel the range is given by 

AK = AL (-rr<K M^)T 


(l«>J 


where Ah is the range of applied loads at a given time, f(— ) is the 


finite width correction factor, (a, b) is the tip stress reduction 

K 

factor, and b is half the stringer spacing. For a fixed value of 'a’, 
(3-3, 

— is a function of the random load parameters AL and r. Thus at a 
dN 

given crack length say a - a^, the growth rate is a random variable. 
The expected value of teh growth rate is given by 



where f (r , AL) are the density functions of the random variables 

r andAL respectively. ^ is the range space forAL and r 

respectively. Equation (11), thus gives the mean crack growth rate at 
any value of crack length under the random loading- This quantity 
expressed in terms of the discrete length units ^ is required in the 
equation for P(k, t) of previous section 
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The mean crack growth rate as given by Equation (11) is a compli- 
cated integral to be solved and does not have a closed form solution. 

Hence, numerical methods have been used to solve the equation. However, 
for a special cases where r atidAl* are stationary Gaussian processes, 
Taylor's series expansion has been used to obtain approximation. Then 
E[a3 at any value of ai is given by the following equation. 

+A \ -f- } /f B 6 "A } H- J 

Mternate Method of Estimation of Parameters 

An alternate method of estimation of parameters is to use the 
fatigue failure data from the same fleet, similar fleet or from tests. 

Such a method requires the following steps. 

(1) The first required step is the solution of equation (7) and 
(8) to obtain P(k, t). This could be left iii the form of 
quadratures . 

(2) The next step needed is the normalization to a realistic maximum 
crack length N (Ah). 

(3) Next step is to estimate the parameters. This has been done 
by using maximxim likelyhood method and by using loss function 
concept in decision statistics. 

This work has been carried out as a part of the project investigation. 
Preliminary results are published in references 11 and 12. These papers 
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A 


xticlude fche consideration of data from a specific fleet supplied by NASA. 

1 ? 

detailed analysis including the model verification will be published 

Applications of the Developed Model 
One of the applications of the developed model is to develop a 
reliability-based, cost-effective design procedure. This method has 
been developed and reported by the investigators in reference 18, 

Some of the significant items and example problems are discussed here. 
Problem Setting 

The problem setting can be best explained by considering an example. 

In this report, the deisgn of a built-up structure such as a sheet- 
stiffener combination is considered. Figure 3 illustrates the stiffened 
panel. The panel is of width w and thickness t. The panel is assTimed 
to be made of a specific material and the particular structure is assumed 
to be a sub-assembly of an aircraft structure. It is also assumed that 
large number of aircraft will be produced as a result of this design. 

Even though the discussed methodology considers a specific material, an 
optimum choice among several candidate materials can be made by following 
a similar procedure and statistical decision theory. External loading 
F consists of a sustained loading Fj^, and a random fatigue loading 
It is assumed that the random fatigue loading has been quantified 
probabilistically. Thus, the total loading F is specified prohahilistically 
For a particular choice of the thickness t, the stringer Spacing 
2b, and the choice of the material, the initial ultimate load carrying 
capacity F^ is known. If the initial microsized flaws or cracks are 
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specified by a probabilistic distributioa, tfea initial load carrying 
capacity is characterised by an appropriate probabilistic distribution 
which depends on the initial flaw size distribution, the material and the 
dependence of the load carrying capacity of the structure on the flaw 
size and other dimensions. 

On the other hand, if it is assumed' that the effect of initial 

flaw size distribution can be described by a crack initiation probability 

distribution^ the load carrying capacity F^ can be expressed as a determinis 

tic quantity if the material properties are also assumed to he deterministic 

The corresponding initial ultimate stress is defined to bec^. Similarly, 

for a given thickness, stress corresponding to external loading is denoted 

by a,. Ifou and (T^ are deterministic, the initial safety margin i.e., 

^ u L 

before fatigue effects are present, is given by the ratio of o' to ct_ 

U Jj * 

As explained earlier, both have. uncertainties and need probabilis- 

tic representation. Then the initial reliability can be considered as a 
safety measure. This can be represented by the probability that cr^/ 
is greater than 1. Due to the presence of fatigue loading, cracks grow 
in size. Crack growth rates and the crack sizes depend on the material 
properties, stress and the number of cycles.' The presence of a crack 
of size reduces the ultimate strength from *^ui‘ reli- 

ability which is defined by the probability that the ratio-^^j^to is 
greater than 1 is also reduced. Consequently, the probability of failure 
which is the probability of the ratio 5: cr^ . The probability 

of failure increases as teh crack lengths increase. The probability 
of failure can be reduced by increasing the initial margin 


of safety or reliability. This of course, increases the weight 
of the structure. Another way of decreasing the probability 
of failure is to inspect the structure at selected times so that the 
cracks can be detected and repaired before they reach their critical 
sizes. In this process, allowable initial margin of safety can be small 
because cracks are not allowed to grow to their critical sises. This 
process however, increases the cost due to inspection. Increasing 
weight also increases the initial cost and the cost of operation. There 
fore, the required design procedure consists of selecting the design 
variables such as the thickness, stiffener spacing, and inspection fre- 
quency during the projected service life so as to minimize the total 
expected cost or weight. The cost and weight can be considered as inter 
changeable functions that can be optimized. Iifeny a time it is easy to 
express the objective function to.be optimized as an equivalent weight 
function. This entire procedure, however, is subjected to the restraint 
that the margin of safety or -reliability does not fail below an accept- 
able limit during the projected life of the structure. 

Therefore, reliability -based fail-safe fatigue design procedure 
consists of selecting specified design variables including inspection 
frequency, subject to constraints, so as to minimize the expected cost 
or weight function while the probability of failure is kept below 
specified limits during the projected life of the structure. In order to 
make the design procedure acceptable to a designer who is not familiar 
with the statistical method.s, the reliability or probability of failure 
can be related to a 'variable’ safety factor or safety margin. 
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Methodology 

The following are the steps that need to be followed in the method- 
ology for the reliability-based fail-safe fatigue design procedure dis- 
cussed in this paper. 

. The first step consists of specifying the design variables and 
constraints. This step identifies the design variables that can 
be selected by the designer to minintize the objective -function 
(weight or cost). 

. The second step is to specify the probabilistic distribution of 
the external loading. This can be a stochastic process, 

. The third step is to formulate the objective function. This 
can be a weight or cost function and is related to the probability 
of failure, the projected life of the structure, the specified 
and selectable design variables, and external loading. 

. The fourth step is to select trial design variables and obtain 
the initial margin of safety or reliability. 

. The next step is to obtain the variation of crack size and crack 
growth probabilities with time. A stochastic model for crack 
growth developed by the authors is used in this report to obtain 
the probabilistic description of crack sizes. This probability 
depends on the material, load description and the number of cycles. 

. from this knowledge of the probability distribution of crack 
sizeg reduction in strength and probability of failure is estimated 
The inspection and repair frequency during the projected design 
life is included in this estimate of the probability of failure, 

13 • 



, The seventh step is to substitute all the information in to the 
cost or weight that was formulated in the third step. This yields 
the cost or weight due to the particular selection of the trial 
design variables. 

. Steps two to seven are repeated with different trial variables to 
minimize the oLiective function by search method. 

. The final design variables are selected subject to restraints such 
as reliability bounds, minimum spacing, etc. 

These are the general steps that are necessary in the design proce- 
dure developed in this report. This needs the description of a stochastic 
model for fatigue crack growth and crack sizes, methods of estimation 
of the probability of failure, methods of including the effects of in- 
spection and repair frequency during the projected design life in the 
probability of failure, and an objective function in terms of cost or 
weight. The stochastic model and the estimation of the parameters of the 
model are already discussed in previous sections. The estimation of the 
probability of failure, reliability improvement due to inspection and 
repair, formulation of the objective function and its minimization are 
discussed in the following sections . 

Probability of Fail,ure 

In this section, method of estimating probability of failure is 
discussed. The improvement in reliability due to inspection, repair and 
consequent renewal and the estimation of this reliability improvement 
are not discussed in this section. These are discussed in the next 
section. 
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The first step in estimation of probability of failure is to 
identify the possible failure modes. In addition to the fatigue failure 
mode, other failure modes such as the sudden over stress or buckling are 
possible. If the event ..of fatigue failure is denoted by the event 

of sudden over stress by E^ and the event of buckling failure by Ej^. 

The probability of failure is given by the union of the all the 
possible events of failure. 

pf ^ o 


Os) 


Probability of occurrence of each of these events depends on the 
strength of the structure to resist ;that particular type of failure and 
the probability of occurrence ot the load that can result in that parti- 
cular type of failure. Because the discussions of the paper are primarily 
restricted to fatigue failure, it will be assumed that only fatigue 
failure are possible. This means that only failure mode possible is 
due to the growth of fatigue cracks and consequent reduction in strength. 

Before discussing the probability of failure under conditions of 
uncertainty, a deterministic design procedure is briefly reviewed here. 
This review is useful in identifying the different probabilistic fatigue 
failure modes. Consider the stiffened panel shorn in Figured. Let it 
be assumed that a central crack is likely to develop in this structure 
due to fatigue. For given w and assumed length between stiffeners 2b, 


the variation of the residual strength a^dwith half the length of the 
29-31 

central crack is shown in Figure 4. The value of the maximum 

external load L is precisely known in deterministic design^ Then, for 
a particular choice of the initial safety margin S, the thickness t and 
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tha corresponding stress cr , critical crack length a can be obtained. 

i C 

These are shown in Figure. 4. As the fatigue cracks initiate and grow, 
failure is not pojssible until the crack attains a length of a^. The 
length of a^ can also be obtained analytically from the following 
formula in the case of a stiffened panel. 



22 

In this equation f(a^/w) is the width correction factor , is the 

23 

tip stress reduction factor , is the fracture toughness of the 
material. 

Because the maximum load L is known precisely in a deterministic: 

case, the stresses due to external load never exceed the residual strength 

for crack lengths a< a . Alternately, it can be stated that probability 

of failure is zero for crack lengths a<a^ and the probability of failure 

is one for a > a . 

— c 

In reality, the external load is not precisely known. The load is 
usually characterized by a random variable. This is the case in which 
reliability based design procedures are needed. In this paper, external 
loading is assumed to be characterized by a stationary stochastic process. 
Even in this case, a value of can be selected in the Figure 4. This 
curve is assumed to be known deterministically. This means that for a 
given width of the panel w and a choice of stiffener spacing 2b, a value 
of critical crack length a is chosen. This value of a corresponds to 
a definite value of on the curve in Figure 4. But, the external load- 

ii 

ing is not known precisely as in the deterministic case. Therefore, the 
value of a and cannot be related to initial safety margin and choice 

C Xj 
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of thickness t. However, the probabilistic description of the external 
loading L is known. As will be shown later, the choice a , cr , and 

C 

thickness t can be related to reliability or probability of failure. 

From a knowledge of the specified bounds on reliability, a and t can 

0 

be chosen. 

Alternately, the following procedure can be used instead of starting 
with a choice, a . A value of (y is selected such that 

5 = = ^ 

whereat (—) is the mean value of external load divided by the choice of 

thickness t and^(^) is the corresponding variance. The quantity is 

constant which is siimilar to safety margin in a deterministic design. 

However, is not arbitrary. The quantities , t and a are related to 

c 

reliability. They can be selected on the basis of the prescribed reli- 
ability bounds. As can be seen in the figure, a selected value of 
corresponds to a value of which corresponds to a value of a^. 

TTnllke the case of deterministic loading, failure may take place 

even for values of crack sizes smaller than a . Such a failure is 

c 

possible because the externally induced stress (L/t) has a probability 
distribution and does not represent the absolute maximum possible 
stress. For values of a<a^, fatigue failure is possible if the externally 
induced stress exceeds the residual strength at any time during 
the service life of the aircraft. This failure is defined as static 
fatigue failure 

In order to simplify the procedure for estimating the reliability, 
the concept of critical crack sise fatigue failure has been introduced. 
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At any time during the service of the aricraft, L«a crack sise at 
a given location can be either greater than or equal to a specified 
length a^ or less than . Then, the following two mutually 
exclusive events can be defined. 

(a) failure occurs when the crack length at the location is 
less than 

(h) failure occurs when the crack length at the location is 
greater than or equal to . 

The event (a) has been denoted by the event (b) will be denoted 

Pfe ■ Pf be Tirrltten as follows. 


^ =7? -p 

•f ^ h 


( 16 ) 


In this equation 


? 


"P (k :>x, V t I «• 


( 17 ) 


In many cases, for practical reasons, a particular value of critical 
crack length can be defined. The structure is considered to have 
failed if the crack length at a given location exceeds this value. 
For example, this absolute critical crack length can be the critical 
length corresponding to the sustained loading F^^ . Then 


•T’. -- .t; -- V o«ii ,t; 

In this case can be called critical siae fatigue failure. 


( 18 ) 


Reliability Improvement Due to Inspectign and Repair 

If no inspections are done during the projected design life, 
the probability of critical crack size fatigue is given by 





fc] ■- J_y 


( 19 ) 


In this equation . is a number that has been defined during the 

normalization of the probability distributions. As pointed out 
earlier probability distributions are normalised to a realistic 
maximum crack length. The reason for normalisation is that cracks 
in reality do not grow to infinite lengths. Then, 


N 

max 


maximum crack length 


M 


The quantity , in equation (19), is the projected design life 
of the sturcture. The probability of critical crack size fatigue 
fialure can be improved due to inspections. This change in 

probability of failure and hence in reliability can be obtained in 
the following way. 

The projected design life is still assumed to be T^^ 

number of hours or cycles. It isasaumed that one inspection is 
done at number of hours or cycles. It is further 
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assumed that this inspection is conducted at = 0»5 T^j . 

At the time of inspection, if cracks of length k(AL, ) k^CiCih) are ob- 
served, the cracks are repairdd. The quantity k^(Al) is the repair 
threshold crack length. It is further assumed that structure is as good 
as new after repair. This means any further crack initiation and growth 
are to be calculated as though the structure is put into service at 
t = T^ and not at t = o. It is also to be noted that only structures 

with k <k<k are repaired because the structures with k(Al)^k CAI) 
r c c ■ 

have failed due to critical sise fatigue failure. It is implicit that 
the cracks of k(i\l) ■«<. k^(Al) are not repaired. 

There is still another qxiantity to be considered. This is the 
probability of detecting a crack by nondestructive inspection techniques 
if a crack exists. In the first step of the derivation, it will be 
assumed that the repair threshold crack length k^(Ab .) is chosen that 
the detection probability is one. Then, the probability of critical 
size fatigue failure in the two intervals can be obtained as follows. 

The probability of failure 1(1) in the first interval corresponding 
to 0<tg,I^ is given by 


By referring to Figure 4, the probability of survival in 0< t< T 


PO) - Jt 


is 1 - P(l) because there is the probability P(l^ that structures fail 
in 0 <t < T . For t^ T , — 


and the probability of repair is given by 

^ = y0 £ 4^ ^ A < i , 
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Then the total probability of critical crack size fatigue failure 

follows; 

p, jp(^>- /0“ 


in 0 <t <2T = T-^ can be written as follows; 

— o D 


,^C = PO) +■ Pg HO + 


where 

and 


F I = P C A , TS J 

p (2-j ^ p 


(‘^^0 


Equation (23 )) for the probability of failure under one inspection is 
obtained by considering the three mutually exclusive and exhaustive 
events F^, and P (1) [see Equation (21 )3. The quantity in the paren- 

thesis of the last term of Equation (23) is the conditional probability 
that the structures will fall in T^< t^2T^ given that they survived 
i:0<^t^T^. This expression for satisfies all the limiting conditions. 
For example, when P^ = 0, P^^ reduces to P(2)9 as expected. 

When Pg^ != 1 and hence P(l) = 0, P^^ becomes zero. Similarly, the 
probability of failure under any number of inspections can be obtained. 

If the crack detection probability due to nondestructive inspection 
techniques is considered, the probability of repair P^ changes. The 
repair is now possible only if a crack of size k^(^l)< k(^l) ■^k^(^l) 
exists and is detected by the HDI capability, with a probability D(k) , 
Here, D(k) is the probability of detecting a crack of size k( 1). 
depends upon the NDI accuracy. A representative function for 

D, is assumed as follows, 
k 


0 

a 4 a. } 


a^ a, 

a„“a. 

2 1 


1.0 



In the illustrative problem, it is assumed that a^ = 0.02” and a^ = 0,3”. 
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Then, the unconditional probability o£ detecting and repairing cracks 
o£ siae k (AI-) ^kCAt) <k (Ab) at T^ is given by 

5 = T.’ m 

r T - , 

Then, of the repairable aircraft given by Plk^^ik-^k^, T^j only Fj^'-'are^ 
repaired and the others are not repaired. "Kow, equations similar to (23) 
can be written with detection probability for cracks included. 

Total Weight Function 

Every optimization problem involves the so-called objective function 

which is a function of the design variables appropriate to the problem 

at hand . The optimum values of the design variables are obtained 

by finding thestationary locations of the objective function subject to 
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the desxgn constraints 

For aircraft structures "weight" is the most crucial consideration 
in design. In the present context, the vreight of the stiffened panel is 
considered to be minimized. The design variables are the thickness of 
the sheet and the width of the stringer spacing. The total "v7eight function" 
comprises of the deterministic weight of the panel and the expected gain 
of weight is given by the product of the probability of failure under 
a given number of inspections and the deterministic weight of the panel. 

The deterministic weight of the panel consists of the weight of the 
sheet and the stringers. Expressed mathematically, the total weight 


function is given by 

W (b/t J = (/+ Pf)[hrtLy + J 

where w = total width of the sheet 
t = thickness of the sheet 
h = breadth of the sheet 
Y*” = density of the sheet material 
N ^ = number of stringers 

St 
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2b = stringer spacing 
^st ~ weight of one stringer 

Equation (26) Is the proper objective function for the minimigation of 
the weight. The effect of increasing the thickness is to reduce to 
expected loss of weight because of the reduction in the probability of 
failure. On the other hand, the deterministic weight is increased by 
increasing the thickness. Thus, a balance has to be found between the 
two. Stringer spacing has the opposite effect on the different weights. 

The minimization is carried out by the search method. The total 
weight function is calculated for a set of thicknesses and stringer 
spacings. It is then plotted versus thickness with stringer spacing a 
the parameter. Then, the lowest weight is selected. The thickness and 
the stringer spacing corresponding to the minimum weight are the optimum 
values if the reliability constraint is satisfied at these values. The 
weight can be expressed in terms of equivalent costs. 

Total Cost function : 

If the problem at hand is the determination of- the optimum number 
of the periodic inspections, then the total weight function may not he 
the proper objective function. Then the total cost function concept 
has to be introduced. The total cost function comprises of the expected 
cost of failure and the deterministic cost of the /periodic inspections. 
The expected cost of failure is given by the product of the probability 
of failure under the given number of inspections and the deterministic 
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cost of structure. The deterministic cost of inspections is proportional 
to the number of inspections. The mathematical expression for the total 
cost function is given as follows s 

G^(j) = + JC^ (27) 

where is the probability of failure under j inspections 

G is the cost of new structure 
s- 

is the cost of one inspection 
J is the number of inspections 

Equation (27) gives the proper objective function because as the ntanber 
of inspections increases j the expected cost of failure decreases while 
the cost of inspections increases. The minimum value of the total cost * 
function is found by the search method. The minimization is subject to 
the reliability constraint,. 

Illustrative Example 

In order to illustrate the developed method, two examples have been 
considered. The first problem is that of a minimum weight design of 
7075-T6 alloy. The problem has been deliberately kept simple for purposes 
of illustration. A more detailed problem is discussed in the Appendix II. 

The design life is supposed to be 15,000 cycles with two periodic 
inspections made during the design life. The reliability is to be 
99.5%. The design variables to be selected are the thickness t and the 
spacing of the stringers 2b. The following data is assumed to be known. 




0.61m 

8.9 

5000 cycles 
175200 N/m 
0;5 


P 

A 

Y = 46.39 Kg/M' 


K, 


IC 


7475/. 23 X 10^ 


2 

<r 5.L 
<r / 


w 


St 


n 


5 X 10 


-13 


0.3048ia 

307 X 10^ (N^/m^ 
0.01 
0.05 Kg 
3 


shape parameter of Weibull distribution for distribution 
for crac'k initiation 

scale parameter of Weifaull distribution. The probabilistic 
description of the loading is given in terms of the mean 
and variance (^2) of the load range and load ratio r. 

As outlined in the preceding sections, the solution procedure is 
carried out. As a first step, the residual strength-critical crack 
length diagrams are obtained for a choice of number of stringers, e.g. 

3, 5, 7j 9, 11, etc. (fig. 6). As the number increases the stringer 
spacing decreases. As one might expect, the rate of growth decreases 
with the nuiiider of stringers. The tip stress reduction factor Cj^(a/b)h 
which is required in the expression for the residual strength is obtained 
from references (29-30) as shorn in Figure 7. 

The variation of the static reliability with residual strength and 
thickness is shotm in Figure 8. For a given loading, in order to 
maintain the same static reliability, the thickness has to increase 
as the design residual strength decreases and vice versa. 



26 


In Figure 9, the relation between the probability of static 

failure, fatigue failure, and total fialure is delineated. 

The total weight functions are calculated in the manner explained 

previously for fixed Rg = 0.9996 and ~ 5j 7, 9, Figure 10 

depicts the minimization curves. From these curves, the minimum. W for 

each curve can be obtained, and then compared with other minima of 

other curves. The overall minum in Figure 10 occurs for a thickness 

of 0 . 106 inches , N .. . = 7 . 

’ strxnger 

Figure 11 represents the minimzation curves for R = 0.9997. As 

expected, the minimum values are now changed, and occur at different 

thicknesses. The minimum now occurs for N . . - 1 and thickness 

strxnger 

t = 0.1044 inches. From Figxjre 12, for R = 0.9998, the overall 

s 

minimum decreases to 3.554 and at N ^ . =7 and t = 0,1052. 

strxnger 

Then the Static reliability R is increased further to R = -.9999 

S' ■ S - . 

the overall minimum is higher than before, i.e. W . = 3.630 and occur 

mxn 

for K ^ . = 7 and thickness t - 0.1052, Figure 13, 

strxnger > & 

Thus comparing all the minima over the various variables, the 

minimum most is W - = 3.554 for R = 0.9997, t = 0.1044 inches and 

mxn s 


N 


= 7. R , the static reliability is equal to (1-P^_,), where 


stringer ' ' "s’ ' ^ ^ '■ fs 

is the static fatigue failure probability. The relation between 

the critical crack fatigue fialure and P ^ the total probability 

of fialure is as follows 


This is plotted in Fig. 9. 

R = 0.9996 does not lead to minimum because (1) the corresponding 
thickness is higher than for R^ = 0.9997. 
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(2) the effect of lower R (0.9996) implies that the structure 

s 

is being designed for a lower residual strength (design stress) from, 
figure 8. This figure is the plot of the following eqtiation not 
quoted in the report. %. 




Ck ^ fP ~'l 


<n.' 




( 




) 




where is residual strength 

is the mean of the load 

2 

is the variance of the load 
t is the thickness. 

The lower design residual stress results in a higher by 
virtue of longer critical crack length (see Fig. 6). The interaction 
of all these results in a higher total weight function for = 0.9996 
It is to be noted that the deterministic weight of the 

panel. This corresponds to an overall reliability of 0.99765 and 
a design residual strength = 15,500 psi. The reliability constraint 
is satisfied since 0,99765 is greater than the reliability bound 


Rjj = 0.995. 

Check on. the Initial Factor of Safety ; 

The mean and standard deviation of the maximum load L , are 

max’ 

obtained* Than, considering different numbers of standard deviations 
above the mean-maximum load L, the initial factors of safety are 
obtained. For example, for one standard deviation above L, the 
initial factor of safety of the optimum design, based on yield 
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strength Is found to be 3.067. ^^hen two and three standard deviations 
are employed, the corresponding factors df safety are 2,60 and 2.32 
respectively. This provides a comparison of the optimum design 
values and the values from equivalent deterministic designs . 

Cost Optimization ; 

To demonstrate cost optimisation, the designed stiffened panel 
is considered. The only variable nov7j is the number of periodic 
inspections or the inspection interval. Since the panel is of a 
given configuration, its weight is fixed. Hence the total cost 
function Equation (34) is the proper objective function to 
be considered in the present context. 

As a first step, the probability of fatigue failure under j 
inspections, j = 0, 1, 2, 3, 4,... etc. is calculated. These 
values are graphically depicted in Figure 13. Corresponding to 
each 6f these numbers of inspections the total cost function 
is calculated from Equation (34), Figure 15. This is repeated 
for various values of the ratio of the cost of one inspection 
to the cost of the structure Cg. When ~ 0.1> the minimum 

occurs for one inspection. Decreasing 0^/Cg to 0.01, 0.005, 0.001 
renders the minimum to occur at two inspections, three inspections and 
four inspections respectively as delineated in Figure 14. 


Alternate Methods of Improving Fatigue Life 
Fatigue Reliability 

The models for fatigue discussed in the preceding sections do not 
apply to cases for which residual stresses are present near fastener 
holes due to a cold working process such as stress coining. The purpose 
of stress coining is _o improve the fatigue life of the structure. A 
simple method of stress coining in aluminum alloy is to expand the 
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fastener hole of the structural member by drawing an oversized 
tnardel hydraulically through the fastener hole, Many. similar 
processes are available for cold working fastener holes. 

Such cold working processes result in a radial flow of the material. 
This results in residual stresses. Residual compressive stresses surround- 
ing the hole provide protection against the fatigue damage by opposing 
the applied tensile stresses, However, as shown in the investigation, 
there is a zone of sustained residual tensile stresses located at a shftt 
distance from the hole. The maximum tensile stress usually occurs at 
the elastic -plastic boundary. Although the tensile stresses are not 
critical in the point of view of fatigue life of the structure, they can 
cause stress corrosion under certain conditions. 

Therefore, the reliability of a stress coined structure needs the 
oDnsideration of both the fatigue improvement and stress corrosion sus- 
tibility. The first step in such a study is to assess the residual 
stresses and stress coining susceptibility in such structures The 
investigations carried out in the project have been published in references 
lA and 15 , 

By-Products From the Project 

As a by-product of the investigations, the following have been 
developed. Ah improved numerical technique was needed in quantitative 
estimation of the parameters of the stochastic model. This has been 
discussed in Appendix III. An application of the technique has been 
done to random vibration problems. The purpose of the application was 
to verify the accuracy of the technique. 

Another by-product is the application of the principles of analysis 
of variance to study the significance of the present methods of grouping 
fatigue failure data. Preliminary work in the field has been discussed 
in Appendix IV, 

30 


Conclusions and Recoiraiiendations 

It has heen demonstrated that an accurate description of fatigue is 
possible by means of stochastic model. A simple model has been developed 
This model can be quantitatively estimated. The model has been applied 
to develop a procedure for a reliabiiity-based cost-effective fail-safe 
design for aircraft structures. In particular, reliability imprcvement 
due to inspection and maintenance has been considered. 

Deterministic design procedures that do not consider the involved 
uncertainties usually result in an over design. This results in an 
increased wight that affects both cost and performance. Furthermore, 
risks involved in a deterministic design are not known. On the other 
hand, the reliability-based design that uses a stochastic model considers 
the uncertainties that are consistent with the model. Risks in a design 
can be assessed. Consistent with the model considered. Such a procedure 
usually results in lower weight than deterministic designs. This results 
in low operating cost and better performance of the aircraft . A very 
costly item in owning and operating an aircraft is the inspection and 
maintenance during the life of the aircraft. As has been demonstrated 
in the project an optimum schemes can he developed by using a 
stochastic model for fatigue and considering the reliability improve- 
ment due to inspection and repair. Methods of including such reliability 
improvement at the design stage has also been discussed. 

The following further investigations are suggested in the point 
of view of the practical application of the developed procedures. 
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1. Development of different types of stochastic models so that 
the user has a choice depending on the particular application. 

It is necessary that all uncertainties be properly Included 
in the model . Different and more accurate methods of quanti- 
tative estimation and verification of the model are needed. 

2. It is also necessary to develop simple optimization techniques to 
include the combination of discrete inspection costs with other 
costs. This is necessary to avoid the difficulty with local 
minimuras.and provide a simple practical procedure, 

3 . The developed procedures should be modified to include multiple 
locations, and multiple cracks. 

4. It appears as though cold working process will be used to 
improve the fatigue life of most existing and future metal 
aircraft. Probabilistic model for failure of such structures 
that includes both the life improvement and the stress corrosion 
susceptibility has not yet been developed. Such models 

are essential to fully take advantage of the cost and weight 
savings potential offered by the cold working processes. 

5. In the point of view of increasing fuel costs, present levels 
of performance can be maintained only by using a material that 
has a higher strength to weight ratio than that offered by 
present aircraft structural materials. Advanced composites 
have such a potential. Mechanical behavior and failure modes 
of these advanced composites are different from that of metals. 
Instead of developing a deterministic design procedure and 

then modifying the procedures to develop probabilistic procedures. 
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reliability'»based design, procedure should he developed from 
the very beginning. By such a process the weight saving 
potential of advanced composites can be explored completely. 
This needs modification of the project to adopt to failure 
modes of composites. 

6, Development of more accurate cumulative damage estimation 

techniques are essential for both metal and composite aircraft. 
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Appendix I 

in this appendix, a method of estimation of the static fatigue 

failure has been discussed. This fatigue failure is possible when 

the external loading exceeds the residual strength of the structure and 

t t 

the crack sise a is less than a^. By defining a quantity S in the 
follox^ing way 


s - Jctl. 


(T\. = 




The probability of static fatigue failure can be defined as the proba- 
bility of s being leas than or equal to 1. Alternately, reliability 
against the static fatigue failure can be defined as R 

=pzs^i2 = 


(a) 


This probability can be evaluated from the following integral if 


the marginal probability density functions given. 

= J' J" l(2/s) dZ dS (3j 


In this equation, f and g are the marginal probability density functions 

of or and C7T respectively, Z is an auxiliary variable and R„ is the range 

space of Z. The integral given in Equation (3) is difficult to evaluate. 

Instead of evaluating the integral of Equation (3), the following 

alternate procedure can be adopted to evaluate the static- fatigue 
29-33 

reliability R * The generalized Ghebychev inequality is employed to 

determine the reliability R • For any shape of density function h(s), 

s 

the probability that the random variable s lies within a range (d -B") 

< s ^ (d is given by the following inequality 

a- 

p/rd-S) e S S (ct-f-S^O ^ l - 'Js 
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in this equation eZ."3 denotes the expectation operat ioiij, 2S' is the 
width of the strip and d is any particular value of s. The lower 
limit of s namely, is unity, i.e.^= d - 1. Substituting these 

limits in Equation (4) , 2 -n 

p[M6^ pi 

Now, recognizing that 

E(S) - S, the mean value of S, and 

9 ^9 

E(S ) = C7g + S 

the equation(5) reduces to the following form after using Equation (2): 


/?c 5. / 






29 32 

For R to be a maximum it is necessary that * 


aa" ^ 


— 

3 


^ ad^ 

From the first of Equation (7), 


^ 0 




+ <rs/{s-\) 


Prom the second of Equation (7) and (8) 

2 

which is negative for all S andCT^. 


(C) 

C^) 

(S) 

(^) 


Substituting for d from Equation (8) i^i Equation (6) it follows that 


(s) 




Os^+{S-0 


(to) 
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Appendix XI 


Numerical Example: 

The proble^i is to design a stiffened panel subjected to a given 
random loading. The panel can have a central crack extending through 
che thickness. Also, the panel will be subjected to periodic main- 
tenance inspections with attendant repairs of the crack when possible. 
Thus, the design variables involved can be categorized as follows; 

(1) Material parameters 

(2) Geometrical parameters and 

(3) Maintenance parameters, ■ 

The design problem therefore consists of (1) selecting the optimum 

material from a given set of different materials, (2) selecting the 
optimuin stringer spacing and thickness, and (3) selecting the optimimi 
number of periodic inspections. 

The following are assumed to be given and the designer has no.;choice 
in these variables 

SI units 

= 20.0" 0.51m 

^ =15.0 0.38 m 

oi = 7.5 7.5 


^ = 48,000 cycles 
= 6.0 X 10^ cycles 


4 

4.8 X 10 cycles 
6.0 X 10 cycles 


= 99.95% 99,95% 

. • V ’ 

In the above set 45 and n> are the overall dimensions of the panel 

The quantity an dyS characterize the Weibull ^uodel for crack initiation 

- 7 " 5 

of 0 . 005 inches . The design life ^ is to he 6.0 x 10 cycles. The 
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reliability restriant should be 0,9995. The material properties 
are as follows: 

For 7075T6 Aluminum Alloy 
= 68,000 lb. /in 
C - 5 X lO"^^ 
n = 3. ■ 

For 2024-T3 Aluminum Alloy 
83,000 Ib/in^^^ 

G = 3 X lO"^^, n = 3 

A computer piogram has been written to obtain the probability of failure 
for each selected thickness, stringer spacing, material and the number 
of periodic inspections !J during the design life. This information is 
later used in another computer program to obtain the expected cost or 
weight function. The design variables that meet the minimum expected 
cost or weight function subject to reliability constraints are selected. 
The following tables illustrate representative results and the selected 
design variables and the material. 

For the first material, i.e., 2026-T3 the overall minimum occurs for 
6 periodic inspections, 3.3” stringer spacing and sheet thickness of 
0.105”. For 707S-T6, the overall minimum occurs for 6 periodic inspec- 
tions, 3.5” stringer spacing and 0-103" thickness when both minimums 
were compared, 7075-T6 has the lower minimum weight at 6 inspections, 

3.3” stringer spacing and 0.103” sheet thickness. Hence, 7075-T6 
would be the selected material. All the details of the calculations 
will be published in the Fh.D. thesis of Mr. B. Uppaluri and in a 
journal. 
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Appendix XII 

Multiple integration with respect to one independent variable was 

needed in integrating the equations (7) and (8) of the main text to 

obtain t) , Such a technique is also needed in many other 

. . 35 

engineer xng problems. Hunter developed a method of numerical multiple 
integration and called it "the integrating matrix method". He applied 
the technique to forced vibration problem of helicopter rotor blades. In 
Hunter’s method^ the derivation of the integrating matrix consisted of 
dividing the range of integration into PT intervals of equal size and 
H4-1 points. At each of the H+1 points, the values of the integrand were 
obtained and represented in a column matrix|f^. The functional varia- 
tion of the integrand in each interval was represented by an r^^ degree 

polynomial. In order to obtain the values of the integral, an r " 
degree integrating matrix|l^ was constructed b using Kewton's inter- 
polation formula. By multiplying the integration matrix and the 
integrand column matrix, the values of the integral x^ere obtained, Fpr 
multiple integration, the integrand matrix was repeatedly multiplied 
by the integrating matrix |l^. For example, 


Improvement 

The mathematical motivation for the improvement is the fact that 
when an 'r^^' degree polynomial is integrated an + 1)^^' degree 
polynomial is obtained. Thus, the improvement suggested is that the 
degree of the integrating matrix be increased by one after each inte- 
gration is a multiple integral. For example, 



^1 X X 

J j j 





o 0 0 

The Improved method was applied for the following problems : 

(i) Multiple Integration of an algebraic function 0< X'<- 20 

(ii) forced Vibration response of a Canilever beam 

(iii) free vibration of Ganl'lever beams 

The results were compared with the exact solutions. 

In the first example, a constant function f(x) = 1,0 was successively 
integrated four times using a second degree integrating matrix and number 
of divisions W =; 20. The percent error ranged from 200.0 at x - 1 to 
015 at x»20,. The improved method was employed with the same H - 20, 
but with integrating matrices of degree 2, 3, 4, and 5 successively. 

The percent error was sero all through the ra?.ige of integration. 

For the forced vibration problem the span was divided into five 
equal intervals (K = 6) and a second degree integrating matrix is 
snployed four times consecutively. The percent error ranged from 6.4 
at 1/5 span to 0,3 at 5/5 span. The improved technique with the same 
H “ 6 but increasing degree of integrating matrix from two resulted in 
a maximtun percent error of only 0.03. 

for the free vibration problem less than 1 percent error in 
natural frequency and/or less mean square error in mode shape was 
obtained at a lower number of spantirise divisions than in the case when 
the integrating matrix was not altered. Also, the mean square error in 


the modeshape compared to the exact mode shape for any mode was less 
in the improved method than in the method of Hunter, 
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The difference between the two methods decreased as the degree of 
the starting integrating matrix is increased. All the results will be 
published. At present the manuscript is being prepared. 




Appendix IV 


Int r oduc 1 1 pfl ' 

Analysis of variance is a means of deteminiing the homogeneity of 
a large collection of data that have been formed by lumping together 
several small groups of data. The small groups are denoted as "sub 
groups" and the variation between them as "variation between subgroups". 
The name? analysis of variance itself stems from an analysis in which 
the total variation in the entire data is partitioned into component 
parts. These components are used to develop a test statistiec. 


The total variation is expressed by the total corrected sum of 


squares j r . e . 


SS-, - 2: Z 




In this equation 

a is the number of treatments 
X^'is the data point 

is the total sum of data points 
N is the total number of data points, and 
n ^ is the number of data points in li *th treatment. . 

The total variation SS^ can be split up into two components as follows: 

SSj = SS^ + SSjj, (2) 

The term SS^ is variation between subgroups and SS^ is variation within 
subgroups . Then, the follox^ing table is constructed to facilitate the 
analysis of variance. 




Source of 
Varlatioa 

Between sub- 
groups 


Sura, of 
Squares 




Degrees of 
Freedom 


a-1 


Mean Sum 
of SqiJares 


SS^ (a-1) 


F 


SS^/Ca-1) 

SS^/CN^> 


Within sub- 
groups 

Total 


SS 


E 


SS„ 


N-a 

W-1 


SSg/(W~a) 


The value in the last column is compared with the critical F 

value at a given percent of significance and degrees of freedom of (a-1) 

and (W-a) respectively. The data is homogeneous if the F value is 

less than the critical F value . 

If the above analysis of variance indicates that the data is non- 

homogeneouSj then it is desirable to find out which of the subgroups 

form a homogeneous set of data. For this purpose, Duncan's multiple 
37 

range test can be employed. It consists of comparing the modified 
difference between the various means (m^^ - ' with the corresponding 

critical valne The modified means are calculated from the following 

eitpressions 

(m. - m.)’ - (m^ - m^.) ‘ 

2r.r. \/^ 

r . -t- r , / 

3 - 1 / 

where r., r. are the number of replications in each group. The critical 

■ 3. . J . 

37 

values can be calculated from Table II of Duncan's paper. Then all 
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the possible groups are subjected to Duncan's test and those groups 
whose modified mean does not exceed the critical value of R' belong to 

IT 

one homogeneous set of data. 

Application 

The procedure that has been discussed in the preceding paragraph 
is used to analyze the fatigue failure data from a specific fleet of 
aircraft. The objective is to investigate if the fatigue failure data 
from several critical regions can be lumped together. If it is possible 
to lump the data together a small number of probability distributions 
can be used to describe the fatigue failure of the entire structure. 

It is also possible to use the system of lumping to do large number of 
inspections at a few representative locations. 

The particular aircraft under consideration has 92 fatigue critical 
regions. Investigations show that the station group (2 to 15), (33 to 38), 
(41 to 46) and (89 to 92) can be lumped together. Analysis of variance 
tests indicate that these subgroups form a homogeneous set of fatigue 
data. The station groups (1-92) ^ (61-70) and (71-80) cannot be lumped 
together because the test results show that their data varies signifi- 
cantly. These results are quantitatively presented in the following 
table. 
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Group 

Variance 

Result 


Within Groups 

Between Groups 


C 2-15) 

0.9490 

0.7785 

No Significant 
Variation 

(33-38) 

0.2680 

0.5972 

No Significant 
Variation 

(41-46) 

0,5229 

0,4460 

No Significant 
Variation 

(41-46 , 
89-92) 

0.'8026 

0.3890 

No Significant 
Variation 

(89-92) 

0.8457 

0,4224 

No Significant 
Variation 

(61-70) 

0.6846 

3.7367 

Data Varies 
Significantly 

(71-80) 

0.6753 

1.8720 

Data Varies 
Significantly 

(33-38, 
41—46 , 
61-72, 
89-92) 

0.7488 

3.8284 

• 

Data Varies 
Significantly 

1-92 

0.7651 

1.4761 

Data Varies 
Significantly 


Complete details will be published iu a Journal. 
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